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The Faddeev-Yakubovsky equations for the a particle are solved. Accurate results are obtained for 
several modern NN interaction models, which include charge-symmetry breaking effects in the NN 
force, nucleon mass dependences as well as the Coulomb interaction. These models are augmented 
by three-nucleon forces of different types and adjusted to the 3N binding energy. Our results are 
close to the experimental binding energy with a slight overbinding. Thus there is only little room 
left for the contribution of possible 4N interactions to the a particle binding energy. We also discuss 
model dependences of the binding energies and the wave functions. 
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I. INTRODUCTION 


In spite of the tremendously increased computational power of today’s super computers, numerical investigations of 
nuclear bound states are still a challenging problem, even for systems of few nucleons. Investigations promise insights 
into the rich structure of nuclear interactions. To this aim one requires reliable solutions of the dynamical equations. 
In this article we would like to present results for the a particle, which are based on realistic microscopic nuclear 
forces including three-body interactions. 

In recent years forces could be adjusted accurately to the huge amount of available NN low energy scattering data 
0,1,1- The overall agreement of the predictions of these model forces with the data is essentially perfect. As a result 
it has been shown that most of the observables in the low energy regime of the 3N continuum could be predicted 
model independently 0 though the interactions themselves are quite different. On the other side, it is known since 
quite a long time that the 3N binding energies (BE) are quite model dependent and, moreover, are generally smaller 
than the experimental value It is assumed that most of this underbinding is due to three-nucleon forces 

(3NF) and modifications of the NN interaction in the presence of a third nucleon. The latter one is of course also 
part of a three-nucleon force mechanism. 

The nature of these 3NF’s is still not completely understood. It is clear that such forces should already arise 
because of the composite structure of the nucleons, what is partially taken into account by allowing for intermediate 
A-excitation. Other mechanisms of various meson-exchange types will also contribute (for a review see jl^ ). In recent 
years there has been new progress in understandin g th e form of nuclear forces, because of the application of chiral 
perturbation theory (yPT) 0 n m, 0 , 0 , 0 10. From this developments one can expect a more systematic 
understanding of the form of NN and 3N forces in the near future. However, yPT implies a priori unknown constants, 
the low-energy constants, which have to be determined from experimental data. The bound states of few-nucleons 
seem to be an ideal laboratory to determine 3NF parameters, as the BE’s are sensitive to the 3N interaction and they 
are expected to be governed by the low-energy regime of nuclear physics |]0. Therefore, the understanding of nuclear 
bound states is an important contribution to the understanding of the 3NF. 

At present these chiral interactions are not as accurate as the traditional, phenomenological NN forces At 

the order of the chiral expansion parameter considered up to now | 0 | , they do not yet describe the NN phase shifts 
with the same accuracy. Allowing, however, for additional fine tuning a high accuracy description can be achieved [p0[ . 
The aim of this article is to pin down model dependences of predictions for the a particle BE and wave function (WF) 
properties. To insure that differences in the predictions are not due to an inaccurate description of the NN system, 
but are due to the more fundamental differences of the interaction models, we restrict ourselves to the traditional 
models in this paper. The techniques developed, however, will help to appl y a lso the upcoming chiral interactions. 
First investigations, using chiral interactions, have already been undertaken [ 0 , | 0 ]. 

Our approach leads immediately to a basic problem. It has been shown that the 3N interaction cannot be determined 
uniquely and, moreover, that each NN interaction has to be accompanied by a different 3NF [0. For the traditional 
NN interactions, there are no 3NF’s available, which have been derived consistently to them. Therefore, we have to 
rely on 3NF’s, which just take parts of the mechanisms into account, which are expected to contribute to the 3NF. 
These models are, for example, the Tucson-Melbourne (TM) | 0 ] and the Urbana IX (Urb-IX) | 0 | 3NF’s. For the 
different NN interactions these models have been adjusted separately to the experimental BE, as described in 
Section HI A. This scheme is justified for two reasons: 
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• It has been shown that many 3N scattering observables in the low energy regime (below ft! 10 MeV nucleon lab 
energy) scale with the BE. This means that predictions for different model Hamiltonians are equal, when the 
models predict the same BE ||^, An adjustment of the 3NF’s exclude model dependences related to 

this phenomenon. We will see that these effects are also visible for the a particle. 


• In the high energy regime (above ft 100 MeV nucleon lab energy) the predictions for 3N scattering observables 
are sensitive to the 3NF showing that the available 3NF models are quite different ||]. This insures that the 
models applied in this paper cover a wide range of possible 3NF’s. Therefore, our results show the model 
dependences of our current understanding of the a particle, which are related to the structure of the 3NF. 


, i 


Bound states of %ht nuclei have been investigated by several groups using different techniques ^^^^ ^ 
1^ , ^ ^ ^ But much of the work is still restricted to somewhat simplified interactions. Perhaps the most 

advanced calculations covering several nuclei have been performed by the Argonne-Los Alamos collaboration ^ . 
Using the Greens functions Monte Carlo (GFMC) technique, they were able to predict BE’s for the light nuclei up to 
A = 8. However, their work is restricted to the AV18 NN interaction model and the class of Urb-IX 3NF’s (new terms 
not considered here have been added in |2^ ). This leads to the question, whether the other available interactions give 
similar or different results for these nuclei. In this respect the “no-core” shell model approach (NCSM) might be 
more flexible. But the work on 3NF’s has not been finished yet. Therefore, we think that a study of the 4N system 
can provide important new information on the nuclear interactions, if one can investigate a wide range of NN and 3N 
models in this system. 

In this paper we use the Faddeev-Yakubovsky scheme to solve the non-relativistic Schrodinger equation (SE) for 


four nucleons. This has been started already in Refs. 


With this method we are able to get 


reliable results for the BE and the WF of the a particle for several NN and 3N interactions. The calculations are 
restricted to A = 3 and A = 4, but we were able to pin down the dependence on today’s interaction models. The 
highly accurate WF, which result from the calculations, are necessary for the analysis of several on-going or planned 
experiments on the a particle, which might reveal the short-range correlations in nuclei or give insights into the 
charge independence breaking of the nuclear interaction Exact WF’s are also necessary to understand the results 
of parity violating e“ scattering experiments Therefore, we will also give first results of calculations including 
the isospin T = 1 and T = 2 component of the a particle ground state WF. 

Another important issue is a first estimate of the size of a possible 4N interaction. We expect that it should show 
up especially prominent in the a particle, because of its high density. Our calculations give some hints, as to whether 
there is room for an important contribution of the 4N interaction in nuclei given today’s NN and 3N interaction 
models. 

In Section H we briefly review the 4N Faddeev-Yakubovsky formalism. The calculations are based on adjusted 
3NF’s. The adjustment procedure is described in Secti on [HI A . Our results for the BE’s of the a particle based on 
various nuclear force combinations are given in Section HIB and the properties of the obtained WF’s are presented 
Finally we summarize in Section 0 - 


in Section HI C 


II. THE 4N YAKUBOVSKY FORMALISM 

The technical challenge in all investigations of nuclear bound states is the accurate inclusion of all short-range 
correlations in the nuclear WF. Due to these short-range correlations, the partial wave decomposition of nuclear WF’s 
is very slowly converging. This hold especially for the very tightly bound a particle WF. Therefore, a rewriting of 
the Schrodinger equation for the 4N system 


^ + E = (1) 

\ i<i i<j<k j 

according to the formalism of Yakubovsky Q is useful. We take NN pair potentials Vij and 3N potentials Uj/c into 
account. T denotes the kinetic energy operator, H the full 4N Hamiltonian and 'k the 4N WF. We will use Jacobi 
coordinates (see Fig. ^ to represent our WF and dynamical equations. These separate the center of mass motion and, 
at the same time, guarantee a kinetic energy operator independent from angular variables. But these coordinates do 
not include all kinds of pair coordinates at the same time and it is hard to describe the short range correlations of 
pairs in other coordinates than their own relative coordinate. Other coordinates unavoidably lead to strong angular 
dependences or, in other words, to a very slowly converging series of partial waves. On the other hand, the Jacobi 
coordinates include the relative coordinates of some pairs. Gorrelations of those pairs are easily described. The WF 
contains the correlations of all pairs and is hard to expand in Jacobi coordinates. This makes the decomposition of 









3 


lA 




FIG. 1: Definition of the lA and 2A type of Jacobi coordinates. 




FIG. 2: The three parts of a meson exchange 3NF, which differ only by an exchange of the particles. 


the WF in Yakubovsky components (YC) highly advisable. The YC’s single out clusters of the four particles. The 
way they are defined guarantees that they are driven by correlations within these clusters only. Therefore, they are 
efficiently expanded in Jacobi coordinates, which single out the same clusters. 

In the isospin formalism nucleons are identical particles. This implies several symmetry properties, which connect 
the different YC’s and reduce the number of independent coupled equations and YC’s to two. The following set of 
Yakubovsky equations (YE’s) are obtained for the two YC’s ipiA and ^2A 

'4’IA = ' 0 ( 12 ) 3,4 = Gotl 2 P [(1 ~ P34) '4’IA + 4’2a\ + (1 + ^ 0 ^ 12 ) Gq C /23 (2) 

tp2A = '1^(12)34 = Goti2P[{l — P34)lplA + 4’2 a] (3) 

The other YC’s are replaced by transposition operators Pij and combinations P = P13P23 + P12P23 and P = P 13 P 24 , 
acting on the two remaining YC’s. The kinetic energy enters through the free propagator Go = 

interaction by means of the pair t-matrix < 12 . The 3NF’s show up in the interaction term Vt 23 - This defines a part of 
the 3NF in the cluster (123), which is symmetric in the pair (12) and which can be related by an interchange of the three 

particles to two other parts and Y /23 that sum up to the total 3NF of particles 1,2 and 3: V 123 = 1 ^ 12 ^ 3 +^ 1*23 + 1^123 ■ 

( 3 ') 

For 3NF’s based on a meson-exchange picture, ¥^23 describes the interaction induced by a meson interchanged between 
particles 1 and 2 and, on the way, re-scattered by the third particle, as indicated in Fig. |[ 

Applying a combination of transpositions to the set of YC’s, one obtains the WF as 

Tr = [1 _ (1 + P)P34](1 + P)^^A + (1 + P)(l + P)^2A (4) 

The YC’s 01A and '02A are anti-symmetric in the pairs (12) or (12) and (34), respectively This guarantees the 
total anti-symmetry of the WF iF. 

The YC’s are expanded in their “natural” Jacobi coordinates. This means that ipiA is represented in the coordinates 
shown in the top of Fig. |^, because both single out the pair (12) and the cluster (123). ip 2 A singles out both pairs, (12) 
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and (34), and is simplest, when expanded in the coordinates, shown at the bottom of Fig. Each of the coordinates 
involves three relative momenta pi2, Ps and 54 or pi 2 , P34 and q, respectively. The angular dependence is expanded in 
partial waves, leading to three orbital angular momentum quantum numbers for each kind of coordinate: ^ 12 , I3 and 
I4 or ^ 12 , I34 and A. We use jj coupling. Therefore, we couple, as indicated in the figure, the orbital angular momenta 
and corresponding spin quantum numbers to the intermediate quantum numbers ji 2 , I3 and I4 or ji 2 and j 34 , and 
these are coupled to the total angular momentum J and its third component M, using two additional intermediate 
angular momenta js and /: {{ji 2 l 3 )j 3 l 4 )JM or {{ji 2 ^)Ij 34 )JM. For the isospin quantum numbers (see Fig. |^) 
similar coupling schemes to total isospin TMt involve only one intermediate quantum number r: ((ti 2 ^)t ^)TMt 
or {ti 2 t 34 )TMT- 

As we already pointed out, the partial wave decomposition requires a huge number of partial waves, whenever one 
needs to represent correlations of pairs and clusters in coordinates, which don’t single them out. Unfortunately, this 
is still necessary in the intermediate states in Eqs. (||) and (||) (for P ipiA, etc.). Therefore, we still need a trem endous 
number of partial waves to find converged results. However, we are going to show numerically in Section III B that we 
can speed up the convergence greatly by using the Yakubovsky decomposition. For our calculations, we decided to 
truncate the orbital angular momenta, requiring that jij < 6 and h, X < 8. Additionally, we constrain the expansion 
for both kinds of coordinates by another parameter requiring that I 42 + I 3 + I 4 < and U 2 + ^34 + A < . 

We use = 14. Our most sophisticated calculations, including the T = 1 and T = 2 isospin channels, need a 
total number of 4200 partial waves for the first kind of coordinates and 2000 for the second kind. We require 36-40 
mesh points to discretize the magnitudes of each of the momenta pi 2 , P3 and <74 or pi2, P34 and q. This insures that 
we obtain results for the binding energy of the a particle, which are converged within 50 keV. 

Using this partial wave truncation, we find that the discretized integral kernel for the set of Eqs. (i and (||) is 
of the dimension 3 • 10® x 3 • 10®. Clearly this can no longer be treated by standard techniques of numerical linear 
algebra, like the QR-algorithm, and one is forced to use an iterative scheme. A Lanczos type method ^ has 
turned out to be very powerful in the past and also here. Succinctly, for an arbitrary A^-component starting vector for 
the unknown amplitude, one applies the kernel leading to a new vector. This is repeated several times by applying the 
kernel always to the new vectors. That set of vectors is then orthonormalized and the unknown amplitude expanded 
into those elements. Inserting this expansion again into the eigenvalue equation Eqs. i) and (^), one ends up with 
a small set of linear algebraic eigenvalue equations of dimension n, where n counts the number of applications of the 
kernel, n is typically 10-20. The energy eigenvalue E, which is buried as a parameter in the kernel is determined in 
such a manner that the eigenvalue of the kernel is 1 . 

Another challenge is the application of the three-nucleon force. In momentum space and partial-wave decomposed, 
this is a huge matrix of typical dimension 5 x 10"^ x 5 x 10“^ for each total 3N angular momentum and parity. In 
case of the a particle, the 3N subsystem total angular momenta have to be taken into account up to Instead 
of preparing these matrices, we handle the 3N forces differently. They can be naturally broken up into a sequence 
of pseudo two-body forces with a change of Jacobi momenta in between (transpositions). This has been described, 
for the first time in The generalization to the 4N system is described in Appendix This technique is much 
more efficient and even allows one to evaluate the 3NF’s in each new iteration of the kernel - no storage of huge 
intermediate matrices related to 3NF’s is required. 

A typical run on a massively parallel T3E with 128 processors takes 2 hrs to get one eigenvalue and the corresponding 
eigenvector. For our method of parallelization we refer to 


III. RESULTS 


A. Adjustment of 3NF’s 


In this paper we restrict ourselves to the modern realistic NN interactions, which are all fitted to the NN data 
with the same high accuracy and also provide a nn force, which predicts a reasonable nn scattering length. These 
interactions are the AV18 and the CD-Bonn [Q. Additionally we show results for the Nijm I, Nijm II and Nijm 93 
interactions [^, which are not adjusted to the nn scattering length and in case of Nijm 93 give a slightly less accurate 
fit to the NN data. The results for the ®He and ®H BE’s are shown in Table ^ They are based on calculations, 
which take two-body angular momenta up to J 12 = 6 into account and are converged up to 2 keV. The full charge 
dependence of the interaction as well as the n-p mass difference are considered. Also the Coulomb force is included 
exactly as described in 


As well known 
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_Ifll all NN model interactions lead to an underpredicted 3N BE. The underprediction is 

strongly model dependent and ranges from 0.8 MeV to 0.5 MeV for the most modern interactions (see Table Q) though 
their description of the NN data is comparable. For benchmark purposes, we also show results for the expectation value 
of the kinetic energy. These tend to be smaller for the non-local interaction Nijm I and CD-Bonn. This behavior can be 
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interaction 

Eb T 

®He 

Eb T 

AEb 

CD-Bonn 

-8.013 37.43 

-7.288 36.62 

0.725 

AV18 

-7.628 46.76 

-6.917 45.69 

0.711 

Nijm I 

-7.741 40.74 

-7.083 40.01 

0.658 

Nijm II 

-7.659 47.55 

-7.008 46.67 

0.651 

Nijm 93 

-7.668 45.65 

-7.014 44.79 

0.654 

Exp. 

-8.482 — 

-7.718 — 

0.764 


TABLE I: 3N binding energies Eb for different NN interactions compared to the experimental values. Results are shown for 
^Hj^He and their binding energy difference AEb- Additionally, we show the kinetic energies T. All results are given in MeV 


traced back to the softer repulsive core of non-local NN forces. We also show the binding energy difference of the two 
mirror nuclei AEb- One sees that all models underpredict the experimental value. The deviation is somehow larger 
for the Nijmegen interactions, which do not describe the nn scattering length correctly. The additional differences for 
the Nijmegen interactions are, therefore, likely a result of an inadequate description of the NN scattering data. We 
will address the issue of the ^H-^He binding energy difference in Ref. |^; therefore, we do not want to go in details 
here. 

Two possible dynamical ingredients are still missing in our calculations: relativistic effects and 3NF’s. We will 
not address the interesting question of including relativity in few-nucleon dynamics here. Attempts to understand 
this issue can be found in Refs. II, H, II n, 1^ . The results of those calculations are varying. Whereas 
approaches based on held equations, like Bethe-Salpeter or Gross equations generally predict an increased binding 
energy compared to the non-relativistic solution, the calculations based on a relativistic Schrodinger equation predict 
a decreased binding energy. In the latter case the relativistic effects are driven by boost properties, whereas in held 
theoretical approaches additional dynamical effects also occur. The magnitude of the predicted effects is of the order 
of 200 keV. The problem is not yet solved. It has also been observed that relativistic effects and 3NF effects are 
related and cannot be separated in held equation approaches H|. In this paper we neglect all relativistic effects, 
hoping that part of them are included in effective 3NF terms. 

The knowledge regarding 3NF’s is similarly scarce, as for the relativistic effects. It has been shown in Ref. 
that 3NF’s are not dehned independently from the accompanying NN interactions. Two 3N Hamiltonians based on 
two different, but phase equivalent NN interactions, can be augmented by a properly chosen 3N interaction to be 
equivalent in the 3N system. In we formulated the more inclusive statement that one could, in principle, always 
hnd NN interactions, which replace a 3N interaction completely in a 3N Hamiltonian. Ref. H] does not conclude that 
this is always possible. It is clear anyhow that the transformation are complicated and therefore it is not practicable 
to use them to get rid of the 3NF’s. As soon as one includes relativistic features the Poincare algebra inevitably 
enforces 3NF’s |pa, which cannot be transformed away. 

In view of this connection of 3NF models and NN force models, a phenomenological approach to the 3NF is justified: 
given a 3NF model, one adjusts its parameters in conjunction with one NN interaction model to 3N or other nuclear 
data leading to different parameter sets of the 3NF for each NN interaction. 

For the Urb-IX 3NF the parameters have been fixed in conjunction with the AVIS interaction using the ^H BE and 
the nuclear matter density predicted by this combination Q. The TM force originally has not been adjusted in this 
way. Its parameters have been deduced from model assumptions and using ttN scattering data [|3[ |6|, |6|] . It is clear 
that a complete 3NF based on meson-exchange should include not only tt-tt, but also p-ir, p-p and so on exchanges. 
Attempts to include these processes have been done, but conclusive results, fixing the parameter sets, could not be 
obtained |]^ . Therefore, we assume in our study that we can effectively include the effects of heavier mesons in the 
TT-TT exchange TM model by a variation of the ttNN form factor parameter A. It has been observed ||, H| that the 
^H BE is sensitive to this cut-off. The original value A = 5.8m,n. has been fixed by matching the Goldberger Treiman 
discrepancy |^]. However, as has been argued in the form factors are ill-defined, because they strongly influence 
the long-range part of the 3NF. Therefore an adjustment is justified. We emphasize that the aim of this paper is the 
investigation of model dependences due to the different 3NF’s. To this aim we only require 3NF models, which are 
different and which have a sufficiently rich spin-isospin dependence. An adjustment of the 3NF does not spoil these 
requirements. 

We combined in Refs. |^, H] the available NN interactions with the TM 3NF and tuned A to reproduce the ^H or 
^He BE’s. The resulting A values are shown in Table ||. The table also includes results for a modified TM interaction. 
It has been argued in p| | that the long-range/short-range part of the c-term is not consistent with chiral symmetry. 
Dropping it leads to a changed set of parameters, which we refer to as TM’. The parameters of TM and TM’ are 
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interaction 

A 

Ei^U) £;(®He) 

AEb 

CD-Bonn+TM 

4.784 

-8.478 

-7.735 

0.743 

AV18-tTM 

5.156 

-8.478 

-7.733 

0.744 

AV18-tTM’ 

4.756 

-8.448 

-7.706 

0.742 

AV18-tUrb-IX 

— 

-8.484 

-7.739 

0.745 

AV18-tUrb-IX (Pisa) @ 

— 

-8.485 

-7.742 

0.743 

AV18-bUrb-IX (Argonne) |p8| 

— 

-8.47(1) 

— 

— 

Exp. 

— 

-8.482 

-7.718 

0.764 


TABLE II: 3N binding energy resnlts for different combinations of NN and 3N interactions, together with the adjusted form 
factor parameters A in units of m,r. The binding energies for and ^He -B(^He) are shown. For completeness the 

splitting AEb is also displayed. All energies are given in MeV. 


summarized in Table Xl of Appendix 

The fits have been done using less accurate BE calculations not including the isospin T = | component and not 
including the effect of the n-p mass difference. Therefore, the new results for the BE’s, shown in the table, do not 
exactly match the experimental values. The deviations are non-significant for the following study, so we refrain from 
refitting the A’s. We adjusted TM to the BE and TM’ to the ^He BE. The table also shows our results using the 
Urb-IX interaction, as defined in ]24| . 

Table || confirms at the same time a well-known scaling behavior of the Coulomb interaction with the BE of ^He 
The adjusted 3N Hamiltonians predict very similar 3N binding energies and AEb^s. This removes the model 
dependence of AEb found in Table We observe that the model independent prediction for these energy difference 
deviates from the experimental value by about 20 keV. Again, we refer to Ref. for a more detailed discussion 
of this issue. In the same reference, a detailed comparison with hyper-spherical variational calculations is given. In 
Table |, for comparison, we only show the BE’s obtained by the Pisa and Argonne group. We note that the calculation 
by the Pisa group is in full agreement with our results. The small deviation from the Argonne result is not significant 
in view of the comparably large statistical error bar of the GFMC calculation. 

We are now ready to apply the 3N model Hamiltonians, given by the A values in Table |^, to the 4N system. By 
using the models from Table ||, we insure that dependences due to scaling effects, as visible for example in AEb, are 
excluded. Given the very different functional forms of the Urb-IX, TM and TM’ interactions, we can expect to see 
any remaining model dependences in our calculations. 


B. a particle binding energies 


Based on these model Hamiltonians, we solved the YE’s and (g) with no uncontrolled approximation. The 
following results are based on a partial wave decomposition truncated using = 14. It has been verihed that this 
is sufficient to obtain converged BE’s with an accuracy of 50 keV. The binding energies given were found varying the 
energy parameter in Eqs. (^) and (||) until the eigenvalue 1 appears in the spectrum of the set of YE’s. 

Independently, one can check the results with a calculation of the expectation value of the Hamiltonian. We 
emphasize that this is an important feature of our method, which minimizes the possibility of errors in the codes or 
unexpected numerical difficulties. 

For these checks one faces the problem to represent the WF with high accuracy. We already pointed out that 
the WF of the a particle is extremely slowly converging, because there is no set of Jacobi momenta suitable to 
describe the short range correlation in all NN pairs. In Table HI we exemplify the convergence behavior of the WF 


for the AVIS interaction. The normalization and the expectation values of the kinetic energy, potential energy and 
Hamiltonian are shown. The WF’s have been derived from the same set of YC’s, using Eq. (Q). The calculation of 
the WF is based on a partial wave decomposition truncated with = 14. In this way we obtained the WF in 
the two different representations, depicted in Fig. |^. For the expectation values shown in the table, we truncated the 
WF in a second step to the partial waves given by the I'^^m parameter in the first column. It turned out that the 
evaluation of the kinetic energy is difficult, because T amplifies the slowly converging high momentum components of 
the WF. The kinetic-energy expectation values, shown in the fourth and hfth columns of the table, do not converge 
within the chosen partial-wave truncation. However, one can rewrite the kinetic energy using Eq. (^) and the fact 
that the transposition operators commute with the kinetic energy and simply result in a sign change, if applied to a 
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jmax 
^ sum 

< tflT >1^ 

< Tit? >2^ 



T(rnix) 

< T|U|T >1^ 

< T|U|T >2^ 

< H 

< H >2^ 

2 

0.9117 

0.9084 

61.27 

62.14 

91.80 

-110.20 

-110.44 

-18.40 

-18.65 

4 

0.9662 

0.9582 

79.10 

76.11 

96.85 

-117.55 

-118.12 

-20.70 

-21.27 

6 

0.9820 

0.9766 

86.36 

83.49 

97.56 

-120.66 

-120.71 

-23.09 

-23.15 

8 

0.9927 

0.9890 

92.41 

90.09 

97.75 

-121.43 

-121.39 

-23.67 

-23.63 

10 

0.9961 

0.9939 

94.59 

92.93 

97.79 

-121.84 

-121.84 

-24.05 

-24.05 

12 

0.9982 

0.9969 

96.10 

95.04 

97.80 

-121.97 

-121.96 

-24.16 

-24.16 

14 

0.9990 

0.9986 

96.51 

95.70 

97.80 

-122.03 

-122.01 

-24.23 

-24.21 


TABLE III: Convergence of the a particle WF for different truncations of the basis states. The superscripts lA and 2A indicate 
the type of Jacobi coordinates employed. The results are based on a calculation using the AV18 NN interaction and no 3NF. 
See text for details. 


interaction 

Ac H 

T Vnn 

Nijm 93 

-24.53 -24.55 

95.34 -119.89 

Nijm I 

-24.98 -24.99 

84.19 -109.19 

Nijm II 

-24.56 -24.55 

100.31 -124.86 

AV18 

-24.25 -24.23 

97.80 -122.03 

CD-Bonn 

-26.26 -26.23 

77.15 -103.38 

CD-Bonn |3^] 

-26.4(2) — 

— — 

Exp. 

-28.30 — 

— — 


TABLE IV: a particle binding energ y p redictions Ea of several NN potential models compared to the experimental value 
and the “no-core shell model” result |^. The expectation values of the kinetic energy T, the NN interaction Vnn and the 
Hamiltonian operator H are also shown. All energies are given in MeV. 


fully anti-symmetrized WF: 


< >= 12 < '^\T\^ia > -b6 < T'|T|i/;2a > 


(5) 


The right-hand side involves mixed matrix elements with the YC’s. The first term has to be evaluated in the 
lA representation, because ijiiA is given in these coordinates, and the second term in the 2A ones because of the 
coordinates of 11)2A- The results for T based on this equation are shown in the column labeled T(rnix) and show a 
promising convergence behavior. We observe a much faster convergence for the YC’s, which was expected and which 
justifies the YE’s approach to the 4N Schrodinger equation. Based on this experience, we normalize our WF and 
the YC’s using a similar formula for the norm. Consequently, the deviation of directly calculated norms of the WF, 
shown in columns 2 and 3, from one is a measure of the numerical error of our anti-symmetization of the full WF. 

Unfortunately, a similar approach is not possible for the expectation values of the potential. However, the interaction 
does not overemphasize the high-momentum tail and its expectation value is much faster converging. We find a 
reasonable agreement of 0.02 % of the lA and 2A results and convergence of both values to uncertainty of 60 keV. For 
completeness we show the expectation value of the Hamiltonian based on T{mix) and the lA or 2A expectation value 
of V. These values agree within 0.1 %. The expectation values differ from the binding energy result of —24.25 MeV 
by only 20 to 40 keV. This is well within the error of 60 keV, which has to be expected from the convergence behavior 
of V and verifies the accuracy of our results. In the following, we will only present the binding energies, the T(rnix) 
values and the expectation values of H and V based on the lA representation. We consider it more accurate than 
the 2A representation, because the norm is closer to one. 

In Table IV our a particle binding energies are summarized for Hamiltonians based on NN forces only. The results 
are identical to the ones published in Ref. except for AV18, where we present a new calculation, based on a more 
accurate grid and taking T = 1 and T = 2 components into account. Due to the more accurate momentum grid, our 
binding energy changed by 30 keV, well within our estimated numerical error of 50 keV. Therefore, we did not redo 
the calculation for the other interactions. The table also shows a result obtained using the NCSM approach Our 
result agrees with their number within the numerical errors estimated. 

As in the case of the 3N BE’s, the 4N BE’s are also underpredicted by all modern NN force models.The underbinding 
ranges from 2 to 4 MeV, showing that the results are also strongly model dependent. Once again the non-local forces 
predict more binding and, similarity, a reduced kinetic energy. The expectation values of H agree within the numerical 
accuracy of 50 keV with the BE’s E^, which have been directly obtained from the YE’s. 













E(^H) [MeV] 


FIG. 3: Tjon-line: a particle binding-energy predictions FiC^He) dependent on the predictions for the binding energies for 
several realistic interaction models. Predictions of interaction models without (crosses) and with (diamonds) a 3NF are shown. 
The experimental point is marked by a star. The line represents a least square fit to the predictions of models without a 3NF. 


interaction 

Ea 

H 

T 

Vnn 

VsNF 

CD-Bonn-|-TM 

-29.15 

-29.09 

83.92 

-106.16 

-6.854 

AV18-kTM 

-28.84 

-28.81 

111.84 

-132.62 

-8.033 

AV18-kTM’ 

-28.36 

-28.40 

110.14 

-133.36 

-5.178 

AV18-kUrb-IX 

-28.50 

-28.53 

113.21 

-135.81 

-5.929 

AV18-kUrb-IX (Argonne) [2^ 

-28.34(4) 

— 

110.7(7) 

-135.3(7) 

-6.3(1) 

Exp. 

-28.30 

— 

— 

— 

— 


TABLE V: a particle binding-energy predictions Ea for the CD-Bonn and AV18 interactions in conjunction with various 3NF’s, 
compared to the experimental value and the Argonne-Los Alamos result. The expectation values of the kinetic energy T, the 
NN interaction Vnn, the 3NF Vsnf and the Hamiltonian operator H are also shown. All energies are given in MeV. 


In Ref. an fascinating linear correlation of the a particle and BE’s has been observed, known as the 
Tjon-line. Our new results confirm this correlation for the newest NN forces. This is displayed in Fig. |3[ One 
sees that all predictions based on only NN forces are situated on a straight line. However, the experimental point 
slightly deviates from this line hinting to dynamical ingredients beyond the NN interaction and the non-relativistic 
Schrodinger equation. We also observe a strong dependence of this result on the accuracy of the NN force. Omitting 
the electromagnetic part of the AV18 NN interaction leads to 16 keV overbinding for the deuteron. A calculation 
based on this potential resulted in a visible deviation from the Tjon-line. 

In the next step we also include 3NF’s into our Hamiltonian. As discussed above, we adjusted these force in 
conjunction with the different NN interactions. We expect a mnch smaller dependence of the BE’s on the 3N 
Hamiltonians in this case, because we remove in this way model dependences, which are correlated to the 3N BE. As 
one learns from the Tjon-line these are the dominant ones. Our results are given in Table ^ Again we obtained an 
accuracy of the BE’s of 50 keV. The convergence is slower for these calculations. Therefore, we do not find the 
same accuracy for the expectation values as for the BE’s. For these we estimate an error of 100 keV, which is still 
within 0.2 % of the kinetic energy. 

For the NN and 3N forces used, we observe a small overbinding of 60 to 800 keV. These results are also included in 
Fig. I For the TM’ and Urb-IX results we find only small deviations of our results from the Tjon-line. For TM we 
see more deviations. The TM force seems to destroy the correlation between the ^H and a particle BE’s. Though the 
TM’ force and the Urb-IX interaction are quite different, their BE predictions seem to be comparable. Unfortunately 
the adjustment of the 3N force has not been done with the same accuracy for TM’. In view of the very expensive 
calculations necessary to improve the TM’ results and in view of the expected agreement of the TM’ and Urb-IX BE’s, 
we did not re-calculate for TM’, but omit its results in the following argumentation. The average BE for the a particle 
using only a NN interaction (based on the restricted choice shown in Table IV) is -24.9 MeV or 88 % of the a particle 
BE. Based on the TM and Urb-IX results in Table we estimate an average 3NF contribution to the a particle 
binding of 3.9 MeV or 14 % of the experimental BE. From the same results we find an average overbinding of 500 keV 
or 2 % of the BE. The contribution of the 3NF is strongly dependent on the NN interaction due to the adjustment of 
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interaction 

r=or=ir=2 

AV18 

99.992 0.003 0.005 


TABLE VI: Contribution of different total isospin states to the a particle wave function. The values are given in %. 


these forces to the 3N BE. The model dependence of the overbinding is much smaller, but depends on the NN and the 
3N force. One can consider this overbinding as the effect of a missing repulsive 4N force. The average size of this force 
can be expected to be 2 % of the BE in the a particle. Certainly the size of this force will be related to the NN and 3N 
forces used. The approach employed in shows that these 4N forces are related to the 3N Hamiltonians in the same 
way as the 3N forces to the 2N Hamiltonians. We conclude from our results that we have found numerical evidence 
that 4N forces are, indeed, much smaller than 3N forces, at least in conjunction with today’s NN and 3N interactions. 
We do not exclude that new additional 3NF terms could be found, which reduce the necessary contribution of 4N 
forces. The results support the generally accepted assumption that meaningful nuclear-structure calculations can be 
performed utilizing bare NN and 3N interactions in a microscopically self-consistent manner. We expect that 4N 
forces probably show up in heavier nuclei in the same order of magnitude (2 % of the BE). We, therefore, suggest to 
take an error of this size into account, when one discusses BE’s for systems with A > 4, based on present NN and 3N 
forces. 


C. Properties of the a particle WF 


Besides the BE’s, we are also interested in the WF of the 4N system, because it serves as input to several analyses 
of experiments involving the a particle. Most of these calculations are based on plane wave impulse approximation 
(PWIA). These calculations are directly sensitive to the WF. Model dependences of the WF are hints to model 
dependences of these observables. However, because WF’s are not observable themselves, we emphasize that these 
dependences might disappear once the full dynamics are taken into account. 

We start with a contribution of the different isospin states to the WF. Because we made a full, isospin breaking 
calculation for AV18 only, there is only one result shown in Table VI, The results for the 3N system do not depend 
on the interaction used [p5[ . Therefore, we do not expect model dependences here. 

One sees an extremely small contribution of the T = 1 and T = 2 component to the WF. However, it is of interest 
that our T = 1 probability, based on realistic nuclear forces, is larger than the one estimated in There the T = 1 
admixture has been found to be about 7 • 10“^ % and the T = 2 state has not been considered. We found the T = 2 
component nearly twice as large as the T = 1 admixture. Moreover the form of our T = 1 state will also be different 
from the one in p^ . As a consequence the isospin admixture correction to the asymmetry as given in will change. 
A renewed evaluation of that correction, also including the larger T = 2 state has not been carried through, but 
appears interesting in view of ongoing experiments. 

WF properties are also important for comparisons to other calculational schemes for treating the 4N system. Am ong 
the most simple of these properties are the S, P and D-wa.ve probabilities of the WF. Th ese are given in Table VH 
for the models based on the CD-Bonn and AV18 interactions. The values given in Table VH are based on overlaps 


between the YC’s and the WF’s, similar to those for the kinetic energies. These numbers are more accurate than the 
results given in However, the differences are not significant, as they affect only the last digit of the results. 

As expected the orbital S'-state is dominant. The H-state probability is sizable, very similar to results for ^H [^ , 
and the P-state gives only a small contribution. The P-state probability for 3N Hamiltonians, based on CD-Bonn, is 
smaller than those for models based on the AV18. This is related to the smaller tensor force of non-local interactions. 
It sticks out that all 3NF’s lead to an increase of the P-wave probability by a factor of 2. 

This raises the question, whether the considered 3NF really act differently in the 4N system. Because of the scarce 
knowledge on 3NF’s, this issue is very important. It insures that we get insight into possible impacts of 3NF’s in 
general, only if our models cover a wide range of interactions. To verify this issue, we decompose the WF’s into parts 
with different total orbital angular momentum, namely S, P and P-states. Based on these components, we calculate 
the expectation values of the Urb-IX and TM 3NF’s for three different WF’s. One is based on the AV18 interaction 
only, one on the AV18-|-Urb-IX and the last on the AV18-I-TM. Four kinds of matrix elements dominate the total 
expectation value of the 3NF : the diagonal S'-S'-state and D-D-stdXe matrix elements and the overlaps of S'-state with 
P- and P-state. Table VHI shows our results. In the second and fourth columns expectation values for Urb-IX and 
TM are shown for the same WF, based on AV18. One observe a strong disagreement of these matrix elements. The 
diagonal elements for Urb-IX are strongly repulsive. They seem to be driven by the isospin and spin independent, 
phenomenological short-range core of the Urb-IX model. In strong contrast, the S-S matrix element contributes most 
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interaction 

S 

P 

D 

CD-Bonn 

89.06 

0.22 

10.72 

CD-Bonn-|-TM 

89.65 

0.45 

9.90 

AV18 

85.87 

0.35 

13.78 

AV18-tTM 

85.36 

0.77 

13.88 

AV18-tTM’ 

83.58 

0.75 

15.67 

AV18-tUrb-IX 

83.23 

0.75 

16.03 


TABLE VII: S, P and D state probabilities for the “^He wave functions. All probabilities are in %. 


3NF 

Urb-IX 

TM 

WF 

AV18 

AV18-bUrb-IX 

AV18 

AV18-bTM 

S-S 

3.16 

2.74 

-2.34 

-4.09 

S-P/P-S 

-0.96 

-2.10 

-1.22 

-3.56 

S-D/D-S 

-5.44 

-7.46 

2.08 

-0.14 

D-D 

0.59 

0.85 

0.01 

0.06 


TABLE VIIL Contribution of different total orbital angular momenta in the wave functions to the expectation values of the 
Urb-IX and TM 3NF’s. All energies are given in MeV. 


of the attraction in the case of the TM. The attraction of Urb-IX is contributed by the S-D overlap. This is a major 
difference in the action of both models in the 4N system. It insures that we used, indeed, very different 3NF models 
though both are based on the 2tt exchange mechanism. Additionally, we see in columns three and five of the table 
the expectation values based on WF’s for the full Hamiltonian. These expectation values differ sizably from the ones 
based on the AV18 WF. We confirm for both 3NF’s that a perturbative treatment of them is impossible. For the 3N 
system this was already emphasized in Refs. Especially interesting is the S-D overlap of the TM force. 

The AVIS WF result is strongly repulsive, whereas the full calculation leads to a slightly attractive contribution. This 
suggests interesting changes of the 3N configurations in the a particle due to this force. 

Are these changes in the 3N configuration visible in momentum distributions? We start in Fig. ^ with a comparison 
of the nucleon momentum distribution 

D{p) = ^ < 4- J = 0 M = 0| 5(p - 94 ) |4' J = 0 M = 0 > ( 6 ) 

47r 

for WF’s based on different NN interactions. The momentum distributions are angular independent. We only consider 
the T = 0 components here. Therefore, the proton and neutron distributions are equal. Because we include in both 
calculations a 3NF, the WF’s are the result of calculations, which roughly give the same BE’s. This insures that we 
do not find differences, which can be traced back to a higher density of the nucleus. The distributions are equal for 
momenta below p = 1 fm“^ for both WF’s. For momenta between p = 1 fm“^ and p = 2 fm“^ the deviations are 
moderate. Above this momentum the AVIS WF is much bigger. We find a clear difference between CD-Bonn and 
AVIS in this momentum region. 

We do not see similar deviations comparing the momentum distributions for different 3NF’s. This is shown in 
Fig. I The WF’s shown there are based on the same NN interaction, AVIS, but differ in the 3NF used. Again the 
BE’s are comparable and there can be no deviations expected because of density differences. In fact one observes that 
the distributions are nearly equal for all models in the whole momentum range. This indicates a remarkable stability 
of momentum distributions with respect to the 3NF. This is in accord with the same independence of the 3NF choice 
for T, the second moment of the momentum distribution, as shown in Table 0. 

The correlations of two nucleons in nuclei are of great theoretical interest. Defined as the probability that two 
nucleons have a certain distance inside the nucleus, one finds very similar correlations for nuclei with different A 
The correlation is characterized by the strong short-range repulsion of nuclear forces, leading to a small 
probability that two nucleons are close to each other. However, the quantitative results depend on the force model used. 
Eor the 3N system this has been shown in [H, and we find similar results for the 4N system ||5^. These correlations 
are not observable. Therefore, difference in this WF property might not show up in observables. Electron induced 
scattering experiments, which intent to see these correlations also see effects of meson exchange currents (MEC’s) and 
final state interactions (FSI’s). Therefore a complete dynamical description of these processes is necessary 
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FIG. 4: Nucleon momentum distributions in ^He on a logarithmic scale. The distribution functions are based on calculations 
using the AV18+TM (solid line) and CD-Bonn+TM (dashed line) potentials. The functions are normalized to J D{p)dp = 



FIG. 5: Nucleon momentum distributions in ^He on a logarithmic scale. The distribution functions are based on calculations 
using the AV18+TM (solid line), AV18+Urb-IX (dotted line) and AV18+TM’ (dashed line) potentials. The functions are 
normalized to J D{p)dp = 


Nevertheless, we want to show those correlations here. Two-nucleon knock-out experiments are expected to provide 
information on relative momentum distributions (see for instance [@ 0 Ell)- In the PWIA these are sensitive to 
the distribution of relative momenta in the nucleus. Consequently, we show in the following momentum correlations, 
defined as 

rp- 

= ^ <^-J = 0M = 0| S{pi2-p) \S Ms >< S Ms\\-f J = 0 M = 0> (7) 

47r 

They are the probabilities to find a pair of nucleons in a spin state \S Ms > and with a relative momentum p. A 
similar definition in configuration space is given in Ref. , where it has been observed numerically that this function 
has a simple angular dependence, which can be expanded in two Legendre polynomials Pf{p- Cz) for / = 0 and / = 2: 

(p) + c^Ms (p) (p . ) (8) 

It only depends on the angle between the momentum and the quantization axis Cz- In Appendix ^ we give an 
analytical proof of this relation. 

The probabilty for two nucleons to be in a fixed spin state S is given by 


Ns = '^ < 'if J ^ 0 M ^ 0\ S Ms >< S Ms \'i' J = 0 M = 0 > 

Ms 


(9) 


For completeness these values are given in Table IX. In the following we will always normalize the correlations to 
47 r J dp C{p) = 1. The probabilties show the importance of individual channels to the total correlation. 
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interaction 

S = 0 

S = 1 

GD-Bonn 

44.60 

55.40 

GD-Bonn-|-TM 

44.98 

55.02 

AV18 

43.07 

56.93 

AV18-bTM 

42.95 

57.05 

AV18-bTM’ 

42.05 

57.95 

AV18-bUrb-IX 

41.87 

58.13 


TABLE IX: Probabilities Ns to find NN pairs in spin S = 0 and 5" = 1 states in ^He as given in Eq. (^). All probabilities are 
given in %. 



FIG. 6: Spin-averaged NN momentum correlations in '‘He for the AV18+TM (solid line) and the CD-Bonn-|-TM (dashed line) 
interactions. The functions are normalized to J C{p)dp = 


In Figs. ^ and ^ spin independent momentum correlations are shown, which have been obtained by summing over 
all SMs states. Obviously, because of no fixed quantization axis, they are angular independent. The first figure shows 
the momentum correlation for the CD-Bonn-|-TM and AV18-I-TM interactions. Similar to the distribution functions, 
they show discrepancies above p = 1 fm”^^. In contrast, we did not find a similar model dependence for different 3NF 
forces. This is shown in Fig. ^ and suggests that 3NF models to not affect observables, which are considered to be 
sensitive to NN correlations. A search for kinematical regions, where FSI’s and MFC’s are suppressed, might reveal 
these correlations. In this case they should show up for momenta greater than p = 1 fm“^. 

We also show the angular dependence of these momentum correlations. In Figs. ||andH both parts of the correlation, 
as defined in Eq. (||), are displayed for S' = 1 and Ms = 0. The angular dependent part does not depend on the 3NF, 



FIG. 7: Spin-averaged NN momentum correlations in '‘He for the AV18-I-TM (solid line), the AV18-|-Urb-IX (dotted fine) and 
the AV18-I-TM’ (dashed line) interactions. The functions are normalized to f C{p)dp = 
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p [fm '] p [fm-‘] 


FIG. 8: Angular independent (/ = 0) and dependent (/ = 2) parts of the NN correlations in "^He for spin S' = 1 

and its third component Ms = 0, as dehned in the text, compared for different interactions on a logarithmic scale. The 
correlation functions are based on calculations using the AV18+TM (solid lines) and CD-Bonn+TM (dashed lines) potentials. 
The functions are normalized, such that the angular independent part fulfills J dp C{p) = The magnitude \C\ is shown. 
+(—) indicates positive (negative) C'/= 2 . 




P [fm-'] 


FIG. 9: Same as Fig. except that the correlation functions are based on calculations using the AV18+TM (solid lines), 
AV18+Urb-IX (dotted lines) and AV18+TM’ (dashed lines) potentials. 


but for higher momenta, on the NN interaction. Around p = 1 fm ^ the / = 2 part is comparable in size to the / = 0 
part. In this region one can expect a visible angular dependence of the correlation. This is related to the toroidal 


structures found in configuration space correlations in |75 

In recent years a knock-out reaction on '*He with ^ in the final state has received a great deal of attention 
pl| . It has been shown that this reaction might be sensitive to the short-range correlations in nuclei . 
A first experiment has not shown the expected dip in the cross section , which has been tracked back to effects of 
MEC’s and FSI’s. Ongoing experiments probe this reaction in different kinematical configurations, which are expected 
to be more sensitive to the correlations. The cross sections in the PWIA or in the more reliable Generalized Eikonal 
Approximation approach [§3l are connected to the "^He/^H overlap functions 


T{p) = ^<^'J = 0M = 0| 6 {q4 -p) \4>t jtmt >< (pt J = 0 M = 0 > 


( 10 ) 


The momentum of the fourth particle is fixed to p and the state of the other three is projected on the triton state (pt 
with spin jt = ^ and third component TO(. Because of the sum over different orientations of the state, T is angular 
independent. One can show that this is still true, if one fixes rrit 0|. The probability to find a inside the a particle 
is given hy Nt = J dp T{p). For completeness we give our results for Nt in Table The results depend slightly on 
the interaction model, but are of the order of 80 %. Thus, one observes a definite change in the 3N configuration in 
the presence of the fourth nucleon. 
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interaction 

iVt(‘‘He) 

GD-Bonn 

84.46 

CD-Bonn -l-TM 

83.49 

AV18 

82.40 

AV18+TM 

80.84 

AV18-bUrb-IX 

80.33 

AV18+TM’ 

80.54 


TABLE X: Normalization constants Nt of the ®H-p overlap distributions in '^He. Results are given in %. 



1 2 3 4 5 

P [fm-‘] 


FIG. 10: ®H-p momentum distribution T in '^He on a logarithmic scale. The distribution functions are based on calculations 
using the AV18+TM (solid line) and CD-Bonn+TM (dashed line) potentials. The functions are normalized to J T{p)dp = 


Fig. 1^ shows the dependence of T on the NN interaction. The function exhibits a dip structure around p = 2 fm“^. 
The structure is the result of a node in the momentum space s-wave function of the fourth particle relative to the 
other three. This node is a necessary consequence of the short-range repulsion. Parity and angular momenta for 
and ^He guarantee that only the s-wave contributes to T. The figure shows that T, indeed, depends on the NN 
interaction. 

The comparision in Fig. [l| of the results for different 3NF’s show that T does not depend on the 3NF’s. Therefore, 
our results confirm that the measurement of T might be valuable to pin down the correlations of two nucleons due to 
different NN forces (if FSI and MEG effects would be negligible). 



FIG. 11: Same as Fig. except that the distribution functions are based on calculations using the AV18-I-TM (solid line), 
AV18-l-Urb-IX (dotted line) and AV18-I-TM’ (dashed line) potentials. 
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IV. CONCLUSIONS AND OUTLOOK 

We solved the Faddeev-Yakubovsky equations for the bound 4N system in momentum space and obtained converged 
results. Two-nucleon interactions, by themselves, underbind the a particle and leave room for considerable model 
dependences. Taking properly adjusted 3NF’s into account, one can considerably reduce the model dependences of 
the BE’s. The combinations of NN and 3N forces lead, in general, to a small overbinding, suggesting that 4N forces 
are repulsive and much smaller than 3N forces. 

We also investigated model dependences of the WF. For momenta below p =1 fm“^ we do not observe any model 
dependences in the momentum distributions and correlations. For higher momenta, only effects of the NN interaction 
show up, because the 3NF’s do not affect these single nucleon and NN properties. 

In contrast, we found a huge effect of 3N forces on 3N correlations visible in the matrix elements of the 3N force. 
These effects require further visualization in future studies. We also found that that the a particle ground state is 
an extremely pure T = 0 isospin state. The admixtures of T = 1 and T = 2 states are of the order of 0.003 % and 
0.005 %, respectively. This sharpens and questions the result found before in Q. 

These calculations provide a baseline for the analysis of experiments involving the a particle, which require highly 
accurate WF’s and insight into NN-force model dependences. The technical developments presented are also important 
for further studies of nuclear interactions based on yPT. First studies have already been started ||l^, ^ |^. yPT 
allows a systematic derivation of 3NF’s, which are consistent with the NN forces. An investigation of these 3N forces 
requires accurate techniques for solving the 3N and 4N Schrodinger equation, in order to fix the parameters of the 
force and to see their effects. The bound states are an intersting object for these studies, because they are the physical 
quantities very sensitive to 3NF effects and are dominated by the low-energy properties of the nuclear interaction. 
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APPENDIX A: PARTIAL WAVE DECOMPOSITION OF CORRELATION FUNCTIONS 

The spin-dependent correlation functions are angular dependent. In momentum space and for a general nuclear 
A-body bound-state with angular momentum J M, it is defined as 

2 

= J M| 5 {pi2-p) IS* Ms >< S Ms\\'i> J M> (Al) 

M 

The operator S{pi 2 — p) IS* Ms >< S Ms\ acts only on the subsystem of particles 1 and 2, i.e., (12). Therefore 
we choose coordinates, which single out this subsystem and denote the coordinates of the remaining particles by 
aA -2 Ja-2 Ma-2, where we have separated the angular momentum quantum numbers. aA -2 also includes the 
motion of the (12) subsystem relative to the A — 2 spectators. The two-body subsystem is described by the usual 
momentum pi 2 and quantum numbers I 12 , S 12 and ji 2 and the third component mj. Resolving the coupling of the 
angular momemtum of the ( 12 ) subsystem and the spectators to the total angular momentum, one obtains for the 
correlation 

XI H X / ^P 12 Pi 2 J dp'12 p'i 2 ^ 

0:A-2 ll2ll2Sl2s'i2 

JA- 2 Ma -2 jl2jl2™i™j 
{ 312 JA- 2 J, mjMA- 2 M) {j[ 2 JA- 2 J, m'jMA- 2 M) 

< dt JM\pi 2 aA-2 i{h2Sl2)jl2JA-2)JM >< p'12 OiA-2 ((^12®12)jl2'^^-2)>^AI| 4 > JM > 

< Pl 2 {h 2 Sl 2 )jl 2 mj\ 5^{p-pi2) IS” Ms X S Ms\ b'i 2 (^ 12 'Si 2 )jl 2 K > (A2) 

The nuclear bound state WF < ^(2 01 A -2 Wi 2 s'i 2 ) 3 i 2 Ja- 2 )JM\^ JM > is independent of M. We choose M = J in 
these matrix elements and perform the M and Ma-2 summation, using the orthogonality relations for the Clebsch- 
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Gordan coefficients. This leads to 

c ^^ hp ) = E 

CiA-2JA-2 


dpi2 PI 2 J dp[2 p'l2 


h2Sl2jl2 

/' s' 
‘ 12*12 


1 


< Pl 2 (^ 12 Sl 2 )jl 2 TOj| S^{p-pi 2 ) IS* Ms >< S Ms\ bi 2 (^ 12 *i 2 )jl 2 "lj > 


2^12 + 1 

iiCj 

< JJ|pi 2 0:A-2 {{W2Sl2)jl2J A-2)J j ><Pl2 Cy.A-2 {{l'l2Sl2)jl2JA-2)JJ\^ JJ> 
In this form the problem is reduced for arbitrary nuclei to the matrix element 


Ml 2 = 


1 


2ji2 + 1 


< Pl 2 (^ 12 Sl 2 )il 2 Wj | S^{p-pi 2 ) \S Ms >< S Ms\ |Pi 2 (^'l 2 *i 2 )jl 2 Wi >, 


(A3) 


(A4) 


which is diagonal in ji 2 and ruj. 

By inserting the unity operator in states of 3D momentum and resolving the coupling of spins and orbital angular 
momenta, we are able to simplify the expression to 


Ml 




d{pi2-p) d{p[2-p) 


1 


■E 


P 12 P P'i 2 P 2ji2 + 1 

IIVJ 

ih 2 Sji 2 ,mj - Ms Ms) (^'i 2 *S'ji 2 , Wj - Ms Ms) Yi\^^^_M^{p) Yi-^^m^-Msip) 


(A5) 


Using standard techniques, one can recouple the angular momenta to obtain a coupled spherical harmonic yf^,, (pp)- 
It turns out that only p = 0 contributes, which is expected, because fixing the spin only hxes the z-axis. The matrix 
elements depend only on x = p ■ Cz. This dependence can be expanded in Legendre polynomials and one ends up with 


Mi 2 = 


A , A '^(P12 - p) d{p'i2 - p) Y^/_'|S-ji2/_yi2+/'i2-/ / (2^12 + 1)(2?!2 + 1)(2/ + 1) 


S I 12 jl 2 
l'l2 S f 


(S f S, Ms 0) {l,,l[J,00) — Pf{x) 


(A 6 ) 


Because S is restricted to 0 and 1, the order of the Legendre polynomial / can only take the values 0, 1 and 2. Parity 
conservation fixes the phase (—)^i 2 +h 2 = 1 . Therefore, the Clebsch-Gordan coefficient (^ 12 ^ 12/5 00) demands even /’s. 
Because of this, the expansion of the angular dependence contains only two Legendre polynomials: Po{x) and P 2 (x)- 
This proves the form of Eq. (^. From the explicit form of M 12 one also reads off that the Ms dependence is giv en by 
an overall Glebsch-Gordan coefficient. This justifies the fact that we only present results for Ms = 0 in Section III C| . 
We also see that the / = 0 part of C is independent of Mg. Finally, we would like to note that the expressions are 
also valid in configuration space, replacing the momenta by the corresponding distances. 


APPENDIX B: TREATMENT OF THE 3NF EMBEDDED IN THE 4N HILBERT SPACE 


1. TM-like forces 

We consider the 3NF as the successive applications of NN-like interactions, which, however, do not respect parity 
and rotational invariance. Only the full 3NF repects these symmetries |Q. The YE’s || and ^ require the matrix 
elements 


<(I2)3,4|U/|>> 


(BI) 


where we can assume that the state 'L is antisymmetric in the nucleons 123. 
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4(27r)® Fo [m-^] 

a b \m^^] c \m^^] d \m^^] 

TM 

179.7 

1.13 -2.58 1.00 -0.753 

TM’ 

179.7 

-0.87 -2.58 0.00 -0.753 


TABLE XI: Strength constants of the TM and TM’ 3NF’s. The numbers are in units of the nucleon mass mjv = 
938.926 MeV and the tt mass m„ = 139.6 MeV. 



FIG. 12: Symbolic representation of a 3NF, like the TM force, and the definition of the momentum transfers Q and Q ' within 
the two subsystems 


One distinguishes four terms in the TM force, the so called a-, b-, c- and d-term, which are given by their individual 
strength constants. These constants are listed together with Vq in Table XI. 


n 23 = Vo { a Ti-T2 WS, +bT-{-T2 


C fi • To 


(B2) 


where we have separated the isospin operators (Pauli isospin matrices r^) and the spin-orbital operators W. 
The FF’s can be read off from the definition of the TM force in momentum space, as given in Ref. f^. 


Q' +ml Q^ + ml 

Q' +ml Q^ + ml 

W^23 = F{Q'")4f^Q'" 11^3^1 = F{Q^)^-^ 

Q' +rnl Q^ + ml 

= F(Q'^) X Q' 

Q' + ml 

with the momentum transfers Q = ki — and Q ^ = ^ 2 ^ — ^ 2 , as indicated in Fig. |^. The Ci’s are Pauli spin matrices 
and the form factors are chosen to be F{Q^) = 


2 \ O’! ■ Q 


Q 


(B3) 




AMQ2 ■ 

Applied to a state vector -0, all four terms have the form 


^ VF 23 I IF 31 


(B4) 


where we have abbreviated the isospin operators by I. 

By introducing the unit operator in the coordinates, which are natural for the W potentials, we are able to turn 
Eq. (H) into 


< (12)3,4|V''> ~ < (12)3,41(23)1,4'>< (23)1,4'|W'23|(23)1,4" > 
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< (23)1,4"|/|(31)2,4"' >< (31)2,4'"|W 31 1(31)2,4* >< (31)2,4*|4' > . (B5) 

We omit the integrals and sums over momenta and quantum numbers of the intermediate states, in order to simplify 
the expressions and denote by {ij)k, I Jacobi coordinates, which single out the pair ij, the 3-body cluster ijk and the 
spectator 1. 4' originally enters in (12)3,4 coordinates. But because of the antisymmetry of 4' in the (123) subsystem, 
the (31)2,4 coordinates are equivalent in this case. 

Because the W 23 ’s do not respect the symmetries of nuclear interactions, the sum over and "'-states have to 
include other parities or other total angular momenta, depending on whether a/c of 6 /d-terms are considered. The '-, 
"- and "'-sums have also to include unphysical symmetric states of the (31) or (23) subsystems. The matrix elements 
of the coordinate transformations < (23)1,4"|/|(31)2,4'" > are given in Ref. The isospin operator leads to a 

change of the isospin part of the transformation. The new isospin matrix elements have been derived in Ref. for 
the 3N system and are given below for the 4N system for completeness. 

The matrix elements for the different IT’s are summarized below. For the a-term one finds 


< (31)2,4|W3“i|(31)2,4'> 

5{qi-q'i) S{p2-P2) 


9494 


P 2 P 2 


<5/47' <5/2/;(5i4j; Sj , p5 mm ' 1.1 


27/^6 2 2 S31 

i 1 531 2 _ 

(p3i iJ/;,^(p3i,P3i)-P 31 Hi^i{P3i,P3i)) 


^ 31 1 

^31 '^31 J3I 
S 3 I ^31 1 


(B 6 ) 


The operator has no isospin dependence; therefore, it is diagonal in isospin space. The momentum dependence is 
given in terms of the function H, which is a combination of Legendre polynomials of the second kind Qi and their 
derivatives Q). 


1 — TTl^ 

HiiP.P') = ^ - QiiBA)) + QKBa) 


with and Bj, = 


The c-term looks very similar as the a-term and follows, if one replaces H in Eq. (B 6 ) by 


~ ? 7 ? 

Hi{p,p') = -^ {QiiBmJ - Qi{Ba)) - 


ml 


PP 


2 {pp'f 


A2 Q[{Ba). 


(B7) 


(B 8 ) 


For our convenience, we used the abbreviation fc = 2fc -|- 1 in these expressions. The notation of the different quantum 
numbers is an obvious generalization of the notation in Fig. |^. 

Because the W 23 - and Jksi-operators are equivalent up to a renumbering of the particles, the matrix elements are 
equal up to a phase factor 

< (23)1,4|VF2“3''|(23)1,4'> = (-)bi+33i-Ht3i-HiA+4i+*k+i < (31)2,dlWg^f |(31)2,4'> . (B9) 

As we already mentioned, we replace the simple isospin transformation matrix element by a combination of the 
transformation and the isospin operator 


< 


(23)1,4|fi • t2|(31)2,4' >= (“)*^^\/^234i 


11 /' 

2 2 ‘'31 

1 1 i 

2 2 


t23 2 


(BIO) 


The b- and d-terms are a slightly more complicated, because the NN-like potentials are now vector operators. The 
matrix elements of the spherical components + iW^), and W~^ = ^ — iW^) 

decompose into a Clebsch-Gordan coefficient and a reduced matrix element 


< ( 31 ) 2 , 4111 / 3 " 1(31)2,4' >= (J'1J,M'^M) < (31)2,4||1F3i||( 31)2,4' > . (Bll) 


The scalar product in spherical coordinates reads 


11^23 • fG3i = ^(-("IT^" 11 / 37 - 


(B12) 
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Because there is no dependence on the third component of the total angular momentum, neither in the transfor¬ 
mation matrix elements nor in the incoming state, we can analytically perform the M” and /i-summations 


M"ii V J' 


J"-J' 


(B13) 


and recover the conservation of the total angular momentum. The NN-like potentials effectively requires only the 
application of the reduced matrix elements and the additional factor ■ The intermediate states are also 

M independent. 

The generalization of the formulas given in Ref. to the four-nucleon system yields 
< (31)2,4||1T3\||(31)2,4>' 

H4:H4 P 2 P 2 


fs 1 h 
J h J' 

27r 

31^31 ^ 


1 jsi Jai 

h j's h 


2 2 ^31 

1 s' i 

*31 2 




/ I HhApsipY) 

[ S31 S33 I31 J 


and 


: 1 1 

-407rV6 ^ ^ 

^31 ■Sai J 31 
f 6 a 2 


P^l>3l' 


X _ 

\/(2a)!(26)! [ hi I 


(a I hi, 00) {b I hi, 00) 


(B14) 


< ( 23 ) 1 , 4111 ^ 2 ^ 311 ( 23 ) 1,4 >' 




PiPi 


j'a 1 h 

J h J' 


1 J 23 323 
h 33 33 


1 1 323 S 23 


-H i 2407rv/6 (-)^33 y^(-)x ^ | 1 l 1 

X K ^ 



'YpHj{p23P23) X] 


rr,' h 

P23 P 23 


‘23 '®23 323 

h3 S23 J23 
h a 2 


> H12AP23P23) 


1 1 X 


2 2 '^23 

^ ^ 523 


tt-j-fc —2 


^2 Y (2a)!(26)! 1 ^23 ^23 ^ 


2 2 

(a [ 4 , 00 ) iblh 3 , 00 ) 


(B15) 


The momentum dependent functions H and H are given in Eqs. (B7) and (B8). 

Again, there is a simple phase relation between W 2^3 and W^i 

< (23)1,4|1T23'^|(23)1,4'> = (_)bi+s3i+t3i+i3i+4i+‘k < (31)2,411^3^1(31)2,4'> 

The isospin matrix element of the d-term differs from the one for the a-, b- and c-term, given in Eq. (BIO). It reads 

< (23)1,4|t3 • (ti X t3)|(31)2, 4 '> = STT'SMTM!^^TT'24i {—)‘^'^^Jh^i 


(B16) 


y^(-)3A+i 


A i 1 

5 1 ^23 


h t23 


i 1 A 
hi 5 1 


. 


(B17) 
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2. Urbana type forces 

The functional form of the Urbana 3NF is much simpler. One usually expresses the Urbana interaction in terms of 
commutator and anticommutator parts. This reads 

^123 = ^271- [ {-^23j-^ 3l} {^2 • fs, "O • Tl} 

+ ^ [-^23, -^3l] [ t 2 ■ T3, fa • fl] ] 

+c/o r2(r23) r2(r3i) 

The force is explicitly dehned in terms of NN interactions 

Xij is derived from the tt exchange NN force. Therefore, it has a spin-spin part Ui ■ aj and a tensor part 

Sij =3 CJi ‘ Ti (7 j * Tj (Tj • (7 j. 

The radial dependence is given as 


YAr) = 
TAr) = 


m-j^r 


(l - 


1 -f 


+ 




m-j^r 


(l - e—■)’ 


The parameters A 2 ,r, Uq and c for the Urb-IX are given in Ref. [ p4| . 

In Q it is shown that the application of the Urbana force can be rewritten as 

< (12)3,4|V''> = 2^2^ < (12)3,41(23)1,4'>< (23)1,4'|X23|(23)1,4" > 

< (23)1,4"|/-|(31)2,4'" >< (31)2,4'"|X3i 1(31)2,4* >< (31)2,4*|V' > 
+C/o < (12)3,41(23)1,4' >< (23)1,4'|T2(r23)|(23)l,4" > 

< (23)1,4"|(31)2,4'" >< (31)2,4'"|r2(r3i)|(31)2,4* >< (31)2,4*|V’ > • 

The isospin operators are very similar to the ones encountered in the TM force 


(B18) 

(B19) 

(B20) 


(B21) 


(B22) 


I = 2 {fi • f2 - - fa • fl X f2) 
= 2(fi • f 2 + ^ fa • fl X fz). 


It is easy to combine Eqs. ( BIOD and (B17) to find their matrix elements. 

The matrix elements of the NN-like interactions X^i and T‘^{r^i) read in momentum space 


<(31)2,4|X3i|(31)2,4'> 
5{qA-q'A) S(j)2-P2) 


qAq'i 


P 2 P 2 




33133 


Yhiipsi^p'si) ^hii' 


31 


(-3 -I- 4s3i) -f fi^^v (P31,P3i) ^S3i1 ^l3ll'3-,331 


and 


< (31)2,4|r2(r3i)|(31)2,4'> 

^ 6s3^s'^^ fel,p'ai)- 

^4^4 P 2 P 2 

Here the tensor operator can be expressed in simple rational functions of the quantum numbers 


^hll'33331 = ^31 - J31 


^31 = isi — 1 
^31 = J31 
^31 = J31 + 1 


_2 0 6 

2J31 + 1 2j3i + 1 

0 _ 2 0 

G V' 331 (j31 + 1 ) _ (y .731+2 

” 2J31 + 1 ^ ^ 2j31 + l J 


(B23) 


(B24) 


(B25) 


(B26) 


^31 = isi - 1 ^31 = 731 ^31 = 731 + 1 
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We numerically perform the Fourier transformations 

2 

YhAP3i,P3i) = - dr A ii^,{p3ir) Y^{r) ii^Ap'siA 

^ Jo 

/ 2 

7i3iiJ,i(P3i,P3i) = - J dr ii^,{p3ir) T^{r) 

2 

Ti3i{P3i,p'3i) = - dr r^ ji^,{p3ir) T^{r) ii^,{pAr) (B27) 

Jo 

with the usual spherical Bessel functions iiix). 

Because these NN-like interactions are all symmetric with respect to an interchange of the subsystem particles, the 
matrix elements for the X 23 and T‘^{r 23 ) equal those for X 31 and T^(r 3 i), respectively 
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